Relativistic Hartree-Fock model for axially deformed nuclei 
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Axially deformed relativistic Hartree-Fock (RHF) model with density-dependent meson-nucleon 
couplings is established in this work, in which the integrodifferential Dirac equations are solved by 
expanding the Dirac spinor on the spherical Dirac Woods-Saxon base. Using the RHF Lagrangians 
PKOi (i = 1,2,3), the reliability of the method has been illustrated by taking the light ?"Ne, 
midheavy °°Fe, and heavy Pb isotopes as examples. As a preliminary application, the systematic 
study of ?"Ne shows that PKO1 and PKO3, which contain the z-pseudovector (x-PV) coupling, 
improve the description of the binding energy of ?"Ne, as compared to PKO2 and the selected RMF 
Lagrangian. Moreover, it is found that the tensor force components carried by the 7-PV coupling 
can present substantial effects in determining the shape evolution of nucleus. 
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I. INTRODUCTIONS 


In past decades, worldwide developments of the 
radioactive-ion-beam (RIB) facilities and advanced 
detectors [1-5] have largely enriched the field of nuclear 
physics, which has been extended from the traditional 
stable nuclei to the ones far from the stability line in 
nuclear chart, namely exotic nuclei [6-10]. Meanwhile, 
lots of novel nuclear phenomena have been observed when 
approaching the drip lines, such as the quenching of 
traditional magic shells and the emergences of new ones 
[11-19], the dilute matter distributions—halo phenomena 
[20-22], etc. Such interesting novel phenomena not 
only largely promote the development of nuclear physics 
with plentiful new opportunities, but also challenge our 
understanding on nuclear systems from both theoretical 
and experimental sides. 

On the other hand, it is well known that most nuclei 
in the nuclear chart are deformed, except a few of 
them which locate nearby magic numbers. In the early 
1950s, many evidences, such as the relationship between 
the nuclear quadruple moments and shell structure [23- 
27] , the rotational-like spectra [28, 29], etc., indicated 
that nuclei can have shape away from spherical. It is 
worthwhile to mention that some consequences of nuclear 
deformation were even discussed by Bohr and Kalckar in 
1937 [30]. Recently, intensive attentions were paid on the 
evolution of nuclear shape because of the notable progress 
of the laser spectroscopy at RIB facilities [31, 32]. Being 
coupled with the deformation effects, more and more 
novel nuclear phenomena were discovered during the 
exploration of the boundary of the nuclear chart, such as 
the island of inversion [33-36], shape coexistence [37-39], 
the superdeformed and hyperdeformed configurations 
[40, 41], etc. In fact, intensive efforts have been devoted 
to understanding these colorful phenomena, which are 
potentially related to the deformation. 

As one of the typical representatives, the relativistic 


mean field (RMF) theory founded on the meson exchange 
4 diagram of nuclear force [42], also referred to as 
a the covariant density functional theory (CDFT) in 
4s recent years, provides an efficient and predictive tool 
4; in exploring the structure properties of nuclei which 
4 spread over almost the whole nuclear chart [43-49]. 
4a With a renormalized relativistic mean field, i.e., the 
so so-call o-w model [50], the RMF theory presents an 
s: appropriate saturation mechanism of nuclear matter. 
s2 Afterwards, in order to solve the problem of large 
ss incompressibility, intensive attempts have been devoted 
s4 to the modeling of nuclear in-medium effects [51], either 
ss via the nonlinear self-couplings of meson fields [52— 
s 54] or the density dependencies of the meson-nucleon 
coupling strengths [55-58], and the accuracy and model 
reliability in describing various nuclear phenomena are 
s» largely improved with the proposed effective nuclear 
so interactions. With the covariant representation of strong 
e& Scalar and vector potentials of the order of several 
e hundred MeV, the RMF theory provides simple and 
63 efficient descriptions of nuclear structure properties, 
such as the natural interpretation of strong spin-orbit 
couplings [59, 60] and the origin of pseudospin symmetry 
s [61, 62]. 

67 Extending to the nuclei with deformation away from 
e the spherical symmetry, large efforts were devoted to 
solving the partial differential Dirac equations of the 
zx nucleon in the relativistic framework, for instance by 
expanding the Dirac spinors on the analytic harmonic 
oscillator (HO) [63-65] or numerical Dirac Woods-Saxon 
(DWS) [66, 67] bases. For the former, the Dirac 
7 spinor and the relativistic mean field that contains only 
75 the Hartree terms of the meson-nucleon couplings are 
expanded in terms of fermionic and bosonic harmonic 
oscillators, respectively, which works well in exploring 
the structure properties of the axially deformed nuclei 
[44]. However, when extrapolating to the regions far from 
stability, namely exotic nuclei, it met serious difficulties 
in providing appropriate asymptotic behaviors of the 
e wave functions because of the limit of the HO potential, 
s; particularly for halo nuclei. Although such defects may 
s be overcome by considering a large mount of oscillator 
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shells [68], it strongly increases the numerical cost that 
in fact violates the simplicity of the HO base. Another 
alterative recipe is to employ the local-scaling point 
transformation to modify the unphysical asymptotic 
properties, namely the transformed HO base [69-71]. 
With the fast development of computational technology, 
it becomes possible to expand numerically the wave 
functions on the complete set made of the solutions 
of the Dirac equation with a Dirac Woods-Saxon 
potential [72], namely the DWS base [66]. Compared 
to the HO potential, the Woods-Saxon potential [73] 
would decrease smoothly to zero at large distance, 
which is essential to provide appropriate asymptotic 
behaviors of the density distributions, particularly for 
the exotic nuclei.  Practically, the extension of the 
RMF theory—relativistic continuum Hartree-Bogoliubov 
(RCHB) theory [74]—has achieved great successes in 
describing not only spherical exotic nuclei [47], but also 
the axially deformed ones with the help of the DWS 
base [75-77]. As one of the typical examples, a novel 
mechanism—the decoupling between the deformations of 
the core and the halo in ^^Mg-——has been predicted by the 
axially deformed RCHB theory [75]. 


Notice that the Fock terms, the inseparable part of the 
meson exchange diagram of nuclear force, are dropped 
in the RMF theory for simplicity. Thus, due to the 
limitation of the Hartree approach, the important degrees 
of freedom associated with the m and p-tensor couplings 
cannot be taken into account by the RMF approach. 
In particular the important ingredient—the tensor force, 
that plays significant role in nuclear shell evolutions 
(78, 79], symmetry energy [80], and excitations [81-85]— 
is also missing in the RMF scheme. In the past more 
than ten years, comparable accuracy as the conventional 
RMF model in describing nuclear structure properties 
has been achieved by the density dependent relativistic 
Hartree-Fock(DDRHF) theory [67, 86] with the proposed 
relativistic Hartree-Fock (RHF) Lagrangians PKO(i = 
1,2,3) [78, 86] and PKA1 [87]. Meanwhile, with the 
inclusion of the Fock terms, significant improvements are 
also obtained on the self-consistent description of nuclear 
shell evolutions [78, 79, 88, 89], tensor force [80, 90, 91], 
spin and isospin excitations [92-96], effective masses [86], 
symmetry energy [97-99], new magicity [89, 100, 101], 
and novel phenomena [101-103]. Thus, it would be quite 
valuable to investigate the effects of the Fock terms, 
particularly the 7- and p-tensor couplings, in describing 
structure properties of deformed nuclei. 


In fact, there already exists some attempt to extend 
the RHF theory to describe the deformed nuclei. In 
2011, the relativistic Hartree-Fock-Bogoliubov (RHFB) 
model for axially deformed nuclei was established by 
expanding the Dirac spinors and RHF mean field on the 
deformed fermionic and bosonic HO bases, respectively 
[104]. However, due to the numerical complexity induced 
by the Fock terms, the calculations are quite time 
consuming, especially for the rearrangement terms due 
to the density dependencies of meson-nucleon coupling 
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strengths, and the application is limited only for light 
deformed nuclei. On the other hand, it is not quite 
expectable in exploring the light exotic nuclei due to the 
limits of the HO base. In this work, as inspired by the 
successes achieved by the DDRHF theory, the axially 
deformed relativistic Hartree-Fock model is developed 
by expanding the Dirac spinors on the numerical DWS 
base, and as a preliminary application the role played by 
the m-pseudovector coupling in nuclear shape evolution 
is clarified by taking ??Ne as an example. Notice that 
the pairing correlations are still restricted in the BCS 
scheme in this work, which works well in describing the 
nuclei nearby the stability line while less reliable for the 
unstable ones, particularly in exploring the halo structure 
in exotic nuclei [105]. Moreover, for the nuclei with 
super/hyperdeformation, the validity of the expansion on 
the spherical DWS basis should be tested more carefully, 
which is not the focus of the current work. 

The paper is organized as follows. In Sec. II, we give 
the general formalism of the axially deformed RHF model 
based on the DWS base. In Sec. III, the results and 
discussions are presented, including the space truncation, 
the convergence check, and the description of ??Ne by 
PKOi, in which the role of z-pseudovector coupling is 
analyzed. At the end, a summary is given in Sec. IV. 


II. GENERAL FORMALISM 


In this section, we will briefly recall the general 
formalism of the relativistic Hartree-Fock (RHF) theory, 
and the RHF energy functional of an axially deformed 
nucleus with the spherical Dirac Woods-Saxon (DWS) 
base. To give a complete impression to the readers, some 
details related to the density-dependent meson-nucleon 
coupling strengths, the pairing correlations, the DWS 
potentials, etc., will be also introduced. Meanwhile, 
the numerical difficulties in prior RHF calculations of 
deformed nuclei will be stressed to provide the readers 
an overall understanding of the status. 


A. RHF Energy Functional 


Based on the meson-exchange diagram of nuclear force, 
the Lagrangian density for nuclear systems, namely the 
starting point of the theory, can be constructed by 
considering the degrees of freedom associated with the 
Dirac field—nucleon (Y) and meson fields—the isoscalar 
ones, including the o meson (c) and w meson (w"), 
and the isovector ones including the p meson (p^) and 
m meson (7), and the photon field (A^). Among the 
selected degrees of freedom, the isoscalar mesons are 
introduced to take the strong attractions and repulsions 
between nucleons into account, and the isovector ones 
for the isospin-related aspects of nuclear force, and the 
photon field for the electromagnetic interactions between 
protons [51]. Therefore, the Lagrangian density for 
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nuclear systems can be expressed as 
L —V +t Llst Lot c A-T.Z, (1) 


where ES Lagrangians of the free fields .Z; (0 = 


v9, 0, uw", p", T, and A") read as 
Zu — (ið, — M) v, (2) 
Lo =+ 50" nüuo — Imo, (3) 
Ly =- 10", + omo (4) 
£o--iR REY + im. (9) 
dea sont OMe — smi m (6) 
24 =- SEP, (7) 


with the field tensors ()"" = Fw” — Ow”, Rew = 
O"p" — Q" pH, and F"" = Q"A" — Q"A". Considering 
the Lorentz scalar (c-S), vector (w-V, p-V, and A-V), 
tensor (p-T), and pseudovector (z-PV) couplings, the 
Lagrangian density r that describes the interaction 
between nucleon and mesons (photon) reads as 


— T3 


7 a 1 
fı = 5 = GeO — Gul Wy ony TF * Pu — e9" Ay 


fo v 
+ ay 7"? p" 


I atas Fb. (8) 
Myr 

In the Lagrangian densities, M and mg are the masses of 
the nucleon and mesons, and gg (¢ = o, wt, p") and fy 
(9' = p",7) represent the coupling strengths of various 
meson-nucleon couplings. In this paper, we use the 
arrows to denote the isovector quantities and the bold 
types for space vectors. 

Following the standard variational procedure, one can 
derive the Dirac equations for the nucleon, Klein-Gordon 
equations for mesons, and Proca equation for photon 
from the Lagrangian density - as 


( — iy,0" + M + X)v(z) = 0, (9) 

(D m2)o =- god (10) 

(+ m2)u* =+ guide, (11) 
(+m?) — ait ré 0, bore, (2) 
(4+ mà)g — Mene (13) 

0, FY! Sel typ, (14) 


where the square box O = 9,0". In the Dirac equation 
(9), the self-energy € can be obtained following the 
variational principle. 

From the Lagrangian density (1), one can also obtain 
the Hamiltonian via the Legendre transformation. After 
neglecting the time-component of the four-momentum 


carried by the mesons and photon, and substituting the 
relevant equations, the Hamiltonian can be derived as 


H -T HM, (15) 
$ 


where the kinetic energy (T) and potential energy (V4) 
terms read as 


T= [ dis) iy: Y + Mu) (16) 
Vo =; J dadzx' )(x)b(x' )U s Dsv (x )v(z). (17) 


In the potential energy terms Vj, ¢ represents various 
two-body interaction channels, namely the o-scalar 
(a-S), w-vector (w-V), p-vector (p-V), p-tensor (p-T), 
p-vector-tensor (p-VT), m-pseudo-vector (7-PV), and 
photon-vector (A-V) couplings. The interaction vertex 
D; (z, 2’) are of the following form: 


I'5.s —— fle(x)gele (18a) 
Tw-v enu). (Guo) > (18b) 
Dv =(9p uT), $ (IP) a (18c) 
1 z vis 
lr =p (forve), i (foc "POLIS (18d) 
Doyr Szaz (foo Fk) (gen), 
1 
+ (gov), ggg fon 708), (18e) 
=i n 
Pr-Pv = (aT), - (frys WO Pa (18f) 
Ds =F 0-7]. - 79]. (18g) 


After neglecting the retardation effects, namely ignoring 
the time component of the four-momentum carried by 
the mesons and photon, the propagators D, (z, x’) in the 
potential terms Vy read as 


1 e me|e-a| 1 1 


== = 1 
rs A (19) 


je- r'|’ An |æ — x'|: 

To get the energy functional, we restrict ourselves on 
the level of the mean field approach. Such that the 
nucleon field operator 7 in the Hamiltonian (16) can be 
quantized as 


ei te Ci, (20) 


)- iG) 


where the annihilation operators c; are defined by the 
positive energy solutions of the Dirac equation (9), e; 
is the single-particle energy (c; > 0), and w;(x) is the 
Dirac spinor of state 7. It should be noticed that in 
quantizing the nucleon spinor % [i.e., in Eq. (20)], the 
contributions from negative states are neglected, namely 
the no-sea approximation, in consistency with the mean 
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field approach that is often utilized in studying the 
ground state properties of nuclei [43, 44, 47, 106]. With 
the no-sea approximation, the energy functional E can 
be deduced from the expectation of Hamiltonian with 
respect to the Hartree-Fock ground state |®o) as 


E = (®o| H |®o) = (o| T |®o) + $ ^ (Go| V5 |o) , (21) 
ó 


(22) 


where A is the mass number of the nucleus, and |0) 
is the vacuum state. In the two-body interaction part 
Va, the expectation leads to two types of contributions, 
i.e., the direct Hartree and exchange Fock terms. If one 
considers only the Hartree terms, it leads to the so-called 
relativistic mean field (RMF) theory. If both the Hartree 
and Fock terms are taken into account, it corresponds to 
the relativistic Hartree-Fock (RHF) theory. 

With the obtained energy functional, the variational 
operation, 


sle DE J da) (æ)ypi(æ)| =0, (23) 
leads to the integrodifferential Dirac equation as 
[temm QU 


where e; is the single-particle energy of the state i, and 
the single-particle Hamiltonian h contains three parts, 
namely the kinetic terms hk, and the mean potential 
terms including the local contributions h? and the ones 
from the Fock terms AP: 


[o.: p - 49 M] 6(z — a^), (25) 
h? = [Er(£)ys + No(w) + Y°Es(a)] 6(m — a"), (26) 


ia Yalx,x') Yp(x,x') 
Xalz,x') Xp(x,x') 


na — 


(27) 


In the above expressions, h? contains the local scalar 
self-energy Xs, the time-component of the vector one Xo 
and the tensor one X, and h? contains the nonlocal 
mean fields XG, Xp, Yg, and Yr contributed by the 
Fock terms. In considering the density dependencies 
in the meson-nucleon coupling strengths, namely taking 
go (¢ = o, w", p") and fy (9! = p",7) as functions of 
nucleon density py = v»y9w, the variation (23) may lead 
to an additional contribution to the self-energy Xo, i.e., 
the rearrangement terms Up [58, 67, 86, 87]. 


B. Numerical difficulties in RHF descriptions of 
axially deformed nuclei 


In this work, we restrict ourselves with the axial 
symmetry, i.e., focusing on the axially deformed nuclei. 
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Standing on the level of the mean field approach, the 
nucleons are considered as point-like particles moving in 
an axially symmetric mean field. Consistent with the 
axial symmetry and reflection symmetry with respect 
to the z = 0 plane, the complete set of good quantum 
numbers is denoted as (vnm), where v is the principle 
quantum number, m represents the parity, and m is the 
projection of total angular momentum on the symmetric 
z axis. In the following, for simplicity, we use the index 
i to denote the quantum number set (vz,m S 0). The 
Dirac spinor of the nucleon in Eq. (20) is therefore of the 
following form with cylinder coordinate (0, z, v): 


dist z)e n ide 


1 | forlo, zjeme 
urm\Q, Z, RS 5 28 
FIOI UOS igir (o, 2)e' 01/2 PS 
ig (0, z) +12 


In practice it is not quite straightforward to solve the 
integrodifferential equation (24), particularly for the 
axially deformed nuclei. 

Within the RMF approach, the Dirac equation (24) 
is reduced as a partial differential equation for the 
axially deformed nuclei. It can be solved by expanding 
the Dirac spinor on the deformed oscillator base, in 
which the local mean fields were also calculated by 
expansion in a deformed harmonic oscillator base [65]. 
In Ref. [104], a similar attempt has been performed 
to solve the integrodifferential Dirac equation for axially 
deformed nuclei, and both local and nonlocal mean fields 
were calculated via the expansion on axially deformed 
oscillator bases. However, the calculations of the nonlocal 
mean fields are too time consuming, particularly in 
dealing with the rearrangement terms generated by the 
density dependence of the coupling strengths [104], which 
largely limits the extensive study in heavy nuclei. 

In cylindrical coordinate space (0,2,y), the meson 
propagators D, can be decomposed as 


oo 


1 ' ^ [^ kdk 
am ip(—v") MEL 
Ds ~ Agr D oe [ ag(k) 


uoo 


x J (ko) J (ko'ete, — (29) 
where ag(k) = (k? + m2)” ^. Similar decomposition can 
be obtained for the photon D4 if set mg = 0 [107]. As an 
alterative choice, one can solve the Dirac equation (24) by 
expanding the Dirac spinor on some selected base, such 
as the oscillator base [65, 104]. Then, the RHF mean 
fields can be calculated directly in cylinder coordinate 
space with the obtained wave functions. In principle, 
such procedure shall be valid for the RHF calculation of 
the axially deformed nuclei. However, in the expansion 
(29), the exponential term exp | — ag(k)|z — 2'|] will 
induce some singularity, i.e., rather sharp peaks appear 
at z = 2’. This brings additional algorithm difficulty 
for the integration over the z coordinate. Moreover, the 
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Bessel functions J, result in the oscillating and rather 
slowly decaying integral term in Eq. (29). Numerically, 
it makes the precise calculation of the integration over k 
rather difficult. 

In contrast to the cylindrical space, the propagator in 
the spherical coordinate space (r, V, p) can be expressed 
as [106, 108] 


where Q = (0,0), the index A, is used to denote the 
expansion terms of the Yukawa propagator, and R5, 5, 
contains the modified Bessel functions J and K as 


1 1 
Mos Bq Veg eei mare) Fossa mors), (31) 
A 
1 1 r 
A < (32) 


T Om Dy F 1 ptT 


Comparing the expressions in cylindrical and spherical 
geometries, it is clear that the integration in Eq. (29) 
corresponds to the infinity sum over Ay in Eq. (30) for 
given uy. Apparently, it is not an easy task to overcome 
the oscillating integration over k in Eq. (29). 

Due to the mentioned algorithm difficulties, in this 
work, we consider the expansion of the Dirac spinor in 
a spherical Dirac Woods-Saxon (DWS) base [66, 67]. 
Compared to the oscillator one, the DWS base owns 
the advantages in providing appropriate asymptotical 
behaviors of the wave functions. It is essential for 
the reliable descriptions of the exotic nuclei. On the 
other hand, it will be similar to deal with the energy 
functional (21) as the spherical case [106], and the cutoff 
of the DWS base in expanding the Dirac spinor (28) will 
automatically truncate the infinity sum over A, in the 
propagator’s expansion (30), which avoids the divergence 
naturally. 


C. RHF energy functional for axially deformed 
nuclei with the DWS base 


Restricted with the spherical symmetry, the complete 
set of good quantum numbers contains the principle one 
n, the angular momentum j and its projection m, and the 
parity = (—1)! (l is the orbital angular momentum). 
For convenience, the quantum number & is often used to 
denote the angular momentum j and the parity 7, i.e., 
k = j+1/2 for j = 1-1/2 and k = —-(j + 1/2) for 
|j = l+ 1/2, and 7 = (—-1)^sign(«) [108]. The Dirac 
spinor in the DWS base is of the following form: 


1 f Galt) Quo (9, i) 
Vam (a) Ur b 


where the index a represents the set of quantum numbers 
(nk) (njl), and Qym (also referred as Qu) is the 
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spherical spinor [108]. Here, we use l, and lg to denote 
the orbital angular momenta of the upper and lower 
components of the Dirac spinor Yam, respectively, and 
ly +la = 2j. Notice that if the quantity « denotes (jla), 
—kK corresponds to (jla). Considering both positive and 
negative energy states in the spherical DWS base, the 
expansion of the spinor (28) can be expressed as 


Vysm X E» Ca iWam (a e Worm (X 


where the expansion coefficient Ca, [i represents (vzm)] 
is restricted as a real number, and v4, reads as 


1 Gi Qiem 
l » (35) 
Fin um 


T 


(34) 


Wurm = 5 Ch inim = 


with Gik = Don Ching and Fig = » Cui Pnr: 

In the present work, we concentrate on the RHF 
Lagrangian PKOi (i = 1,2,3) [78, 86]. Notice that 
PKO2 contains the o-S, w-V, p-V, and A-V couplings, 
and in addition PKO1 and PKO3 take the 7-PV coupling 
into account. Thus, the RHF energy functional can be 
expressed as 


E EF" +S" (E? + EP), 
$ 


where EF" = (®o|T|®o), and EP and EJ are the 
Hartree and Fock terms of the two-body potential 
energies E; = (®o| V, |®o) with ó = o-S, w-V, p-V, Tt- 
PV, and A-V. 

With the expansion (34), the kinetic energy functional 
E" can be derived as 


(36) 


dG; 


jn =E f eis.[ d 


[2 


+ “Gin — MFin! 


Gir Fi MGix|\, (37) 


where v? (€ [0, 2]) denotes the occupations of the orbit i. 
For PKOi, only the local self-energies Ng and ‘Np 

remain in h? Eq. (26), and the relevant scalar and 

baryon densities, i.e., ps and pp, can be deduced as 


ps = 25 vi(z)v; 


p» = » = 2% 


where P), is the Legendre een Aa Shall be even 
integer starting from 0 because of the parity conservation, 
and the cutoff of Aq is determined naturally by the 
truncation of the DWS base in the expansion (34). The 
expansion term pò can be derived as 


= > v2 D c0 
x ED + 


Qrr2 


even 


FN 


r) Py, (cos 0), (38a) 


vi (a r)Py,(cosd), (38b) 
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where + in squared brackets corresponds to the baryonic 
and scalar densities, respectively, and for the symbol Z 
see Eq. (A4) for details. 

With the expansion of densities (38), the Hartree terms 
of the energy functional E7 can be expressed as 


EP =) rdrEM, s()pM(r), (40a) 
Be. ys T r?drXo*. v (r)p;" (r), (40b) 


where the expansion term of self-energies are given in 
Eq. (A9). For the other two vector couplings, i.e., the 
p-V and A-V couplings, the expressions can be obtained 
similarly as the w-V ones, and the details can be found 
in Appendix A 2. 

For the contributions from the Fock terms to the 
energy functional (21), namely ES , they can be expressed 
in a unified form as 


ZI dv 5 ( (Gin F. les) 
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where ¢ represents the o-S coupling, the time and space 
components of the vector (w-V, p-V, and A-V) ones, 
and the 7-PV one. Here, we introduce R to denote the 
nonlocal density terms, which appear in the nonlocal self- 
energies Yo, Yr, Xg, and Xp as source terms, 
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For the details of the nonlocal self-energies Ye, Yr, Xe, 
and Xp in Eq. (41), the readers are referred to Appendix 
A3. 


D. Density dependent meson-nucleon couplings 


In the current RHF approach, namely the DDRHF 
theory, the nuclear in-medium effects are introduced 
phenomenologically by taking the coupling strengths as 
the functions of density p, [86, 87], i.e 


=94(Po) fa (£), 
—ge (Oe, 


for @ =0-S,w-V, 
for go! —0p; Js fr; 


(46a) 
(46b) 


go(Po) 
Jo! (Po) 


283 


284 


285 


281 


o 


287 


288 


28! 


o 


290 


291 


where £ = py/po, po is the saturation density, and 


1 + bs (6 + dọ)? 
fe (£) —9$3 ES co(£ dt dg)? 


The parameters a, b, c, d in the expressions (46) are 
determined with the parametrization referring to the 
properties of nuclear matter and the selected finite nuclei 
(78, 86, 87]. 

Under the axial symmetry, the coupling constant gg 
can be expressed in terms of the Legendre polynomial as 


(46c) 


even 


) = V 9? (po) Pa, (cos 0), 
Xp 


9o(Pb (47) 


where A, shall be even integer due to the parity 
conservation and its cutoff is decided by the convergence 


requirement. T'he expansion term jg? is calculated by 


1 
9," (po) =| d 
—1 


Because of the density dependence in gg, one has to take 
the rearrangement term “ip into account to preserve the 
energy-momentum conservation [58]. The variation of 
the coupling constants gg can be expressed as 


=) Pak (cos V) 2. 


ai Ap 


(cos#)Py,(cos)ag(os). (48) 


Ope? 
4 Ca i 


óge( Po) -6Ca ED (49) 


Og" $ 
ap Po 
where the fraction term Ape /0C5,i can be deduced easily 


from the expression of p?" (39), and the term agi” /üp,* 
is derived as 
À pr ( 


E vai 


where 0g /Opy is determined by the density-dependent 
s ; A 
form NU and L = V2L+1. Notice that both gj” 


and 0g j / dpe are calculated numerically, in which the 
integration is evaluated with the Gaussian algorithm. 
Due to the density dependencies in the coupling 
strengths, the rearrangement terms, which appear in the 
self-energy XM, can be simply expressed as 


y (x23 da 4 xb sg 
$ 


where the detailed contributions from the Hartree and 
Fock terms are given in Appendices A2 and A3, 
respectively. 
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E. Eigenvalue equations: Dirac equations with the 
spherical DWS base 


Since the Dirac spinors (28) are expanded on the 
spherical DWS base, the total energy of an axially 
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deformed nucleus is expressed as a functional of the 
coefficient set {Ca}. The variation (23) is therefore 
performed with respect to C, ;, which leads to a series 
of eigenvalue equations as 


Notice that the expansion of the Dirac spinor [see Eq. 
(34)] contains N = (Np + Np) x Km terms, and Np, 
Np and Km will be introduced in the following context. 
Therefore, H^" = (HT is a N x N square matrix, 


aa’ 

and Cc = (Cai) is a column matrix with N terms. 
The eigenvalue, i.e., the single-particle energy £;i, can 
be determined by diagonalizing H7™, as well as the 
eigenvector a 

Similar to the single-particle Hamiltonian in the 
integrodifferential Dirac equation (24), the matrix H™™ 
in Eq. (52) consists of three parts, i.e., the kinetic H*™, 


local HP, and nonlocal H7 terms, 
H =H" + HP + HF, (53) 


where the subscript (7m) is omitted. With the spherical 
DWS base, the kinetic term can be expressed as, 


Akin = ni arf Gal É CET MGy| 
T 
ra R 
+F, | ms MF, Jh (54) 
where a = (nk) and a’ = (n'&), and for HE? it requires 


k = K'. For the local H D that contains the Hartree mean 


aa’ 


fields and rearrangement terms, it can be expressed as 


HE, = f ar DWI (6. Fa) 
Xa 
Gu 
x Es +E EP G ) (55) 
à a! 


where $' represents the vector channels, i.e., the w-V, 
p-V, and A-V couplings. For the nonlocal term, it can 
be expressed with the non-local mean field Ye, Yr, Xe, 
and Xr in the energy functional (41), 


HM / drdr' (Ga Fa), 
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In this work, we focus on the RHF Lagrangians PKOi 
(i = 1,2,3) and ¢ represents the o-S, w-V, p-V, n-PV, 
and A-V couplings. 


F. Pairing correlations: BCS method with Gogny 
force 


For the open-shell nuclei, the pairing correlations play 
an important role in determining the properties of ground 
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state, particularly for the nuclei close to the drip line 
[47, 74, 109-112]. For the deformed nuclei, the pairing 
correlation also palys an essential role in the shape 
evolution of nucleus since the single-particle structure 
will be largely changed with the shape evolution. In this 
work, the pairing correlations are considered with the 
BCS scheme. 

For the even-even nuclei, the BCS ground state can be 
expressed as 


m0 


Il (ui + vic] 


i 


|BCS) = cl) |-), (57) 


where v? (€ [0,1]) represents the occupation probabili- 
ties, u? +v? = 1, and i denotes the time reversal partner 
of the state i. From the variation with respect to vj, 


5 (BCS| (H — AN) |BCS) = (58) 
the gap equations can be derived as 
1 m'»0 As 
A= vip : ; 59 
9 2 " (er — A)? + AZ 9) 


where N is the particle number operator, the chemical 
potential A is introduced for particle number conserva- 
tion, A; is the pairing gap of state i, and the pairing 
interaction matrix element can be written in following 
general form as 
ve = T 


TP = (i VPP Wi) — (si VPP i'i’). (60) 

In order to provide a reliable description of the pairing 
effects, we adopt the Gogny force D1S [113] as the 
pairing force to take advantage of the finite range, i.e., 
the natural convergence with the configuration space 
in evaluating the pairing effects. It is essential for 
the exploration of the shape evolution of the nucleus, 
since the single-particle spectrum is largely changed with 
respect to the deformation. The Gogny-type pairing force 


reads as 


VP (rir) = ye elt) us (Wx + ByP? 
x=1,2 

— H4P' — MyP"P"), (61) 
with the parameter uy, Wy, By, Hy, and My (x = 1,2) as 
the finite range part of the Gogny force, and P? and P7 
are the spin and isospin exchange operators, respectively. 
Consistently with the expansion of the propagator (30), 
the coordinate part of VPP is also expanded in spherical 
coordinate (r, V, p), 
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where the radial part reads as 


2 
—(p3 2 2 H 
Vy,a(r, 7’) =e (r*+ Y nx] 2n; Dena 


In order to abbreviate the expression, we rewrite the 
Gogny pairing force as 


V(r,r') 22x M^ (Ay + DP’) 


x=1,2 


5 Vy, A(T, r’) 
A=0 


x 5 Yi (9, e)Yx, (9. e), (64) 

H 
with A, = W,— H,P* and D, = B, — M,P". For given 
orbits 4 and ï’ in an axially deformed nucleus, the pairing 


interaction matrix element V£? can be derived under the 
spherical DWS base as 


VBP = / drdr! Y (KE ine KE uw) 
YS, YE, KE ilk! 
KK! KK! (r,r’) ik, iK J pt 
where KẸ pẹ and KẸ yp read as 
C; Giy >$ FipFy cl 
KS polr) =E, Kf uwr) ===, (66 
ik, i'R (r) Inr? ” ik, ik (r) Inr? C ( ) 


and the nonlocal parts YC, Y”, XC, and X” have the 
following form: 


ESDA ET O a)? 
XA 
-D4(CMwo)|. (67) 
PE =F P nr) + D) (Opa) (68) 
XA 
Xe, =; Y Ver) (AR Dx) (CHp), (69) 
X, 
XL = ip Vyx(r, r’) | (As T Dy) (c 
"E 
-D4(Cibo)]. 0 


In the above expressions, the sum over A shall fulfill A + 
lu +l! to be even. Notice that in deriving V$’, the L-S 
coupling scheme is used [74] and for simplicity we only 
consider the contribution of the main component J = 0 


(J=L+8). 


G. Dirac Woods-Saxon potentials and deformation 
parameters 


In calculating the nuclear structure under the CDFT, 
the self-consistent iteration starts in general from an 
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initial potential. In this work, we take the local Dirac 
Woods-Saxon potentials [72] as the initial ones, which are 
also used in determining the spherical DWS base. The 
only difference is that in the initial DWS potential we 
need to consider the effects of deformation, i.e., to choose 
appropriate initial deformation parameters 69. With the 
underlying iteration, the calculation will converge to a 
local minimum of the energy functional nearby £o. 

In the Dirac equation, the local self-energies consist of 
the vector and scalar terms, i.e., X9 + y? Xs, which are 
also referred to as X4 = X + Ng for convenience. As 
an initial one, the local self-energy X7 and the nonlocal 
ones in Dirac equation (24) are set to be zero, and the 
Woods-Saxon type local X+ can be expressed as follows: 


1 — ag(N — Z)13/A 


UB(r)=+V +V, Tla 

PC) ?1- exp [(r — RẸ)/a%] ie) 
T3 = = 

U(r) = — Vo EEG A +V, (71b) 


1 + exp [(r — Aa 


where the isospin projection operator 73 is defined with 
the convention that 73 |n) = |n) and 7s |p) = — |p), and 
RE =re A/? represent the empirical radius of neutron 
or proton for the nucleus with mass number A, and 
VZ represents the Coulomb potential between protons 
(ra = —1). The parameters in the Dirac Woods-Saxon 
potential, being used in this work to provide the spherical 
DWS base and initial potential, are given in Table I. 


TABLE I. Parameters of the Dirac Woods-Saxon potential, 
and Vo — —71.28 MeV and ao — 0.4616 [72]. Except the 
dimensionless a,, the other parameters are in fm. 


Qv T+ To, a4 a 
Neutron | 11.1175 1.2334 1.1443 0.615 0.6476 
Proton 8.9698 1.2496 1.1400 0.6124 0.6469 


In the Dirac Woods-Saxon potentials (71), the 
Coulomb potential between protons is evaluated by 
assuming the charge distributed uniformly, 


3 r? ) 
2R.  2R3/ 


V. -az( when r « Re, (72a) 


1 
V. =aZ-, when r > Re, (72b) 
r 


where Re = RP= and a is the fine-structure constant 
that represents the coupling strength of Coulomb 
interaction. 

Different from determining the spherical DWS base, 
one has to take the deformation effects into account 
in the initial DWS potentials (71). If restricted with 
the spherical symmetry, one can divide a nucleus into 
spherical shells with equal density at various radial 
distances. Extending to an axially deformed nucleus, the 
surface with equal density becomes an ellipsoidal one, 
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which can be described as 


R(9,) =Ro|1 + Do eosYa« (9,9), 
Au 


(73) 


where Ro corresponds to the radius of spherical nucleus 
with equal volume. Considering axial symmetry and the 
reflection one with respect to the z = 0 plane, R(V, y) 
can be further reduced as 


R(0; B) =Ro|1 qa ier A cos? 0— ij), 


where 8 describes the deformation of a nucleus deviated 
from the spherical shape, and the positive and negative 
values of 6 correspond to the prolate and oblate shapes, 
respectively. In the DWS potentials (71), the empirical 
radii R shall be replaced by RẸ (0, 8), leading to X'? = 


EP (r, ð; 8), 
IE Vs 3 cos? d — i). (75) 
= rp a(i- =) Gs a , (76) 


where ô = 34/5/1678, and the values of rj? 
in Table I. 

For the initial DWS potentials X¥ (r, 9; 8), similar 
expansion terms as the local self-energies in the Hartree 
energy functional (40) can be determined by 


(74) 


RE (0; B) =Rors 


T3 
Ro: 


are given 


(77) 


T 
x^) =| 37? (r, 9) Py (cos 0) sin Vdo, 
0 
where the integration can be evaluated with the Gauss- 
Legendre algorithm. The initial matrix in the eigenvalue 
equation (52) is then obtained by replacing the expansion 
terms Dit + X in HP, Eq. (55) with X2", since the 
kinetic part HE? only depends on the DWS base and the 
nonlocal part H s is set as zero. 
With the expansion of densities (38), it is rather 
straightforward to deduce the intrinsic multiple moment 
Q5, for a deformed nucleus as 


/ 4 Xa 
Qu = pen f arte (r), 


and with obtained quadruple moment Q2, the quadruple 
deformation f can be evaluated approximately by 


vV5nQ» 
3RZA ’ 


(78) 


p= (79) 


where Ro = 1.241/3, 


III. RESULTS AND DISCUSSIONS 


As the first attempt of the RHF calculations with 
the DWS base for an axially deformed nucleus, we 
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first introduced the space truncation in this section. 
As the convergence check of the space truncations, we 
show the test calculations for the light ?"Ne, midheavy 
56Fe, and heavy Pb isotopes. Further we concentrate 
on the structure properties of ??Ne described by the 
RHF Lagrangians PKO;, as compared to the calculations 
with the RMF Lagrangian DD-ME2 [114] and available 
experimental data. 


A. Space truncations 


Within the DDRHF theory, the calculations for axially 
deformed nuclei in this work are performed by applying 
the expansion of Dirac spinor on a spherical DWS base 
[66], i.e., Eq. (34). Related with that, two additional 
expansions shall be considered carefully, namely the 
decompositions of the propagators (30) in spherical 
coordinate (r,9,«) and the expansions of the coupling 
strengths in terms of Legendre functions (47). Thus, 
one needs to handle the truncations related to (nk) 
in Eq. (34), A, in Eq. (30), and A, in Eq. (47). 
In fact, these three truncations are not independent 
from each other. For a given A, in Eq. (47), the 
propagator's expansion A, is automatically truncated by 
the cutoff on (n) in the spherical DWS base. Eventually, 
two independent truncations remain, namely the cutoffs 
on (n&) in expanding the Dirac spinors (34) and A, 
in expanding the density-dependent coupling strengths 
(47). 

For the Dirac spinor Wyrm with axial symmetry, in 
which the parity m and angular momentum projection 
m are good quantum numbers, the expansions (34) are 
carried with respect to the principle quantum number n 
and & quantity on the spherical DWS base [uos]. For 
completeness, both positive and negative energy states of 
the spherical DWS base have to be taken into account. 
One should not mix this with the no-sea approximation 
adopted in the mean field approach, which corresponds 
to neglecting the Dirac sea in calculating the densities 
or currents. Notice that in expanding the Dirac spinors 
(34)-(35) the sum over the principle quantum number 
n is carried beforehand. Thus, for all the included 
k blocks, the numbers of positive and negative energy 
states can be simply fixed as the same Np and Np values, 
respectively. One just needs to consider enough large Np 
and Np values, and such choice does not bring additional 
costs in the calculations, since the complicated nonlocal 
self-energies are independent on the principle quantum 
numbers v or n (see the Appendix for details). 

In general, the Dirac equation with the spherical DWS 
potentials are solved in a spherical box with the size Rm 
To keep the numerical accuracy, the Np value shall be 
modified with respect to Rm. In most of the cases, it is 
accurate enough to set Rm as 20 fm, and Np and Np 
are chosen as 28 and 10, respectively, and the accuracy 
has been verified in the RHFB calculations [67]. For the 
k cutoff, we use Km to denote the number of « blocks 
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included in the expansion (34) of the Dirac spinor Wyrm 
(28). In total, the number of the expansion terms in Eq. 
(34) is then Km x (Np + Np). 

For the open-shell nuclei, one needs to truncate the 
configuration space in evaluating the pairing effects. 
Practically, such truncation corresponds to the maximum 
values of m and v of each m orbit. In this work, the finite- 
range Gogny force D1S is introduced as the pairing force 
such that the calculations can be smoothly converged 
with respect to the truncation. Thus, one can consider 
v and m values as large as possible. For the former, 
it does not increase the numerical burden in calculating 
the Fock terms which are the dominant costs paid in the 
calculations. However, for the maximum m value Mmax, 
the computation costs can be enlarged by times when 
increasing Mmax- 

In order to unify the cutoffs and reduce the numerical 
complexity, we set the maximum m value as Mmax and 
the number of & blocks considered in expanding the 
Dirac spinor Wurmma, is set as Kym,,,,, which gives the 
maximum absolute & value as kmax = Mmax + Ig, 
1/2. Hence, for an arbitrary Dirac spinor Yyrm, the & 
quantities in the expansion (34) include |x| = m + 1/2, 
m -- 3/2, ---, kmax and sign(k) = m * (—1)l*l, m = 41 
being positive and negative parities, respectively. 

For a given nucleus, the max value is usually 
determined by the convergence requirement of the BCS 
pairing. Practically, both the mmax and kmax (or Km,,,.) 
values are decided by a careful convergence check. For 
instance, in the case of large deformation, some orbits 
containing the components with large « values may 
intrude even across the well-known major shells. Thus, 
one needs to enlarge Mmax and Kym,,,, values, which can 
remarkably increase the computational cost. 

For the density-dependent coupling strengths (47), 
the A, value shall be an even integer starting from 0 
because of the axial symmetry and parity conservation. 
In general, it is accurate enough to have four terms 
Ap = 0,2,4,6 for most of the nuclei. Thus, with fixed 
truncation in expanding the Dirac spinor, the cutoff on 
the A, value in decomposing the propagators (30) is 
automatically determined as Ay = 2kmax + Me 


B. Convergence check 


At first, we take the light nucleus Ne as an example of 
the convergence check. Figure 2 shows the binding energy 
[plot (a)] and the quadruple deformation 8 [plot (b)] as 
referred to the m-cutoff Mmax, in which Kma, is fixed 
as 3. In Fig. 2(a), it is seen that the binding energies 
converge when Mmax > 9/2, as well as the deformation 
parameter 8 shown in Fig. 2(b). On the other hand, it is 
also shown that with different initial deformation fo, the 
calculations converge to the local minimum nearby. One 
may notice that with mmax = 3/2, the calculations with 
Bg = 0 and 0.5 give the same binding energies. For ??Ne 
with spherical symmetry, the last occupied orbits are ds /2 
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which exceed Mmax = 3/2. Thus, the calculations with 
Bo — 0 can not converge to spherical shape. Considering 
the convergence of the binding energy and deformation, it 
is precise enough to choose Mmax = 11/2 and Kma — 3 
for ?°Ne. 
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FIG. 1. (Color online) Binding energy Eg (MeV) (a) and 
quadruple deformation £ (b) of ?°Ne calculated by PKOi and 
DD-ME2 with respect to max, in which Km,,,, is fixed as 
3. The open symbols denote the calculations with the initial 
deformation 69 = 0 and the filled ones for fo = 0.5. 


Here, we only show the convergence check with respect 
tO Mmax- For the expansion of the density-dependent 
coupling strength, A5^* has been checked for the nuclei 
with a wide range of mass numbers, and it is enough to 
choose A5^* = 6 for most of nuclei. For the Kmmax value, 
it is accurate enough to fix Kma, = 3 with appropriate 
Mmax Value that is in general decided by the configuration 
space of the pairing correlations. For economical reasons, 
one can carefully choose an appropriate combination of 
the Mmax and Kma, values. 

We further check the convergence for the medium- 
heavy nuclei, here, taking °°Fe as an example. In Table 
II is shown the binding energy Eg (MeV), the total, 
neutron, and proton quadruple deformations 6, n, and 
B, with respect to the m-cutoff max, in which K7yp,,,. 
is set as 3 and fo 0.3. It is seen that for both 
PKO1 and PKO2, the calculations are converged after 
Mmax > 11/2. Even in the cases with mmax < 11/2, 
the results are rather close to the converged ones. It 
is similar to the calculations of ??Ne with 8o 0.0, 


azo Which converge quickly when Mmax > 5/2. This is due 
40 to the reason that with small deformation or spherical 
41 Shape, the contributions from high « blocks contribute 
4» little to the expansion of the Dirac spinor Yyrm whereas 
43 with large deformation, e.g., the calculations of 7?Ne with 
44 Bo = 0.5 in Fig. 1, the components with large « values 
as may have substantial contributions. 


TABLE II. Binding energy Eg (MeV), quadruple deformation 
B (the total), 8, (neutron), 8, (proton) of °°Fe calculated by 
PKO1 and PKO2 with respect to max, in which Kmmax is 
fixed as 3 and the initial deformation fo = 0.3. Notice that 


ER” = —492.26 MeV [115] and 6°*P = 0.25 [116]. 

PKO1 PKO2 
Mmax Ep B Bn Bp En B Bn Bp 
7/2 | —489.9 0.24 0.24 0.24] —489.3 0.22 0.22 0.22 
9/2 | —489.9 0.24 0.24 0.24] —489.4 0.22 0.23 0.22 
11/2 | —490.0 0.24 0.24 0.24) —489.4 0.23 0.23 0.22 
13/2 | —490.0 0.24 0.24 0.24) —489.4 0.23 0.23 0.22 
15/2 | —490.0 0.24 0.24 0.24) —489.4 0.23 0.23 0.22 


For the heavy nuclei, we take the even Pb isotopes 
from 18?Pb to ?!4Pb as examples for the convergence 
check. Table III shows the binding energies calculated 
by PKOi and DD-ME2, and the experimental (exp) 
values [115] and the results calculated with the spherical 
code (Sph.) are also given as the references. For the 
calculations with deformed (Def.) code, Mmax is set as 
15/2 and K,,,.. = 3. Notice that in both Sph. and 
Def. calculations the pairing correlations are treated 
with the BCS method by taking the Gogny force D1S as 
the pairing force. For all the selected isotopes, the Def. 
calculations with £y = 0 give the spherical shape. It can 
be seen that the deviations between the Def. and Sph. 
calculations are rather tiny (dozens of keV), particularly 
for the doubly magic ?°°Pb, which in fact illustrates 
the accuracy of the Def. code in calculating the heavy 
nuclei. For the open-shell Pb isotopes, the accuracies are 
as perfect as the doubly magic ?99Pb, and the relative 
deviations are usually less than 0.0196. For the matter 
radii, a similar accuracy can be also obtained by the Def. 
code, which is not shown for simplicity. 

As a further example of the convergence check, we 
choose the heavier nucleus ???Rn that has the reported 
experimental quadruple deformation as f 0.1269 
[117]. It is worthwhile to notice that for the heavy 
nuclei the truncations of the expansions (34) and (47) 
shall be tested carefully. Here, we take five terms 
Ap = 0, 2, 4, 6, 8 in expanding the density-dependent 
coupling strengths (47) for ??Rn. Using PKO1 and 
ss DD-ME2, the convergence tests with respect to Mmax 
sx are shown in Table IV, including the binding energy 
s7 Ep, matter and charge radii r and ren, and quadruple 
ss deformations (f, Bn, Bp). As shown in Table IV, both 
s9 RHF (PKO1) and RMF (DD-ME2) calculations show 
5 appropriate convergence with respect to Mmax, which 
sı has illustrated the reliability of the expansion on the 


11 


s» DWS base. Moreover, both RHF and RMF Lagrangians 
533 provide a reasonable description on the bulk properties 
su. of ??? Rn, referring the available data. 


C. Description of "Ne 


535 


s6 After a careful convergence check, we performed 
systematical calculations for ?"Ne, a typically deformed 
nucleus. First, we show the constraint calculation with 
respect to the quadruple deformation f for ?"Ne in Fig. 2, 
where the space truncations are set as Mmax = 11/2 and 
Ka, = 3 and the filled cycles denote the local minima. 


s2 For comparison, we employed the RHF Lagrangians 
saa PKOi and the RMF one DD-ME2 [114]. 


537 


538 


539 


540 


541 


-130 FY _ i _ 20 
he: Minas TA KaT 3] Ne 
-135 Es 
c -140 
> 
© 
= -145 
m 
W -150 
-155 
-160 Lr. 
B 
FIG. 2. (Color online) Binding energy Ep (MeV) as a 


function of deformation 8 for ?°Ne. The results are calculated 
by PKOi and DD-ME2 with mmax = 11/2 and Km = 3: 
The filled circles denote the local minima. 


max 


For ??Ne as shown in Fig. 2, the RHF Lagrangians 
PKO1 and PKOS3, which contain the a-PV coupling, 
present rather similar results. In contrast, the 
RHF Lagrangian PKO2 and RMF one DD-ME2 give 
systematically less bound results. Referring to the 
experimental value of the binding energy of ??Ne, it can 
be seen that PKO1 and PKOS3 show nice agreement, 
and the results given by PKO2 and DD-ME2 are less 
bound and the deviations are 5 MeV. It shall be 
mentioned that there exists some discrepancy on the 
binding energies between our results and the ones in 
Ref. [104] (Figs. 1 and 9 therein) For the RMF 
Lagrangian DD-ME2, our results are coincident with 
the calculations using the code developed from Ref. 
[65], while showing about 1.5 MeV less than the one 
in Ref. [104]. Moreover, for the RHF Lagrangian 
PKO2, the calculations in Ref. [104] (Fig. 1 therein) 
only present nearly converged results for ?"Ne, showing 
similar deviation from our calculations as DD-ME2. The 
s possible reason might be due to the fact that stronger 
sea pairing effects are obtained by the general Bogoliubov 
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TABLE III. Binding energies of Pb isotopes with even neutron numbers form N = 100 to 132, calculated with spherical (Sph.) 
and axially deformed (Def.) RHF codes. The results are calculated with the RHF Lagrangians PKOi and the RMF one 


DD-ME2, in comparison to the experimental data [115]. 


A Bas. PKO1 PKO2 PKO3 DD-ME2 
Sph. Def. Sph. Def. Sph. Def. Sph. Def. 

182 | —1411.65 | —1411.76 —1411.80 | —1411.94 —1412.13 | —1412.61 —1412.64 | —1409.81  —1409.88 
184 | —1432.02 | —1431.28 —1431.31 | —1431.78  —1431.98 | —1432.09 —1432.11 | —1429.22 —1429.31 
186 | —1451.80 | —1450.32 —1450.34 | —1451.18  —1451.38 | —1451.07  —1451.08 | —1448.12 —1448.21 
188 | —1471.07 | —1468.94 —1468.97 | —1470.16  —1470.35 | —1469.63 | —1469.64 | —1466.66 —1466.77 
190 —1489.81 | —1487.19  —1487.22 | —1488.73  —1488.92 | —1487.82 —1487.83 | —1484.91  —1485.01 
192 | —1508.10 | —1505.09  —1505.15 | —1506.91  —1507.08 | —1505.68  —1505.70 | —1502.90  —1503.01 
194 | —1525.89 | —1522.68 —1522.75 | —1524.69  —1524.84 | —1523.23  —1523.26 | —1520.65 —1520.75 
196 | —1543.17 | —1539.97  —1540.05 | —1542.06  —1542.19 | —1540.48 —1540.52 | —1538.18  —1538.27 
198 | —1560.04 | —1556.97  —1557.04 | —1559.01 —1559.11 | —1557.44 —1557.49 | —1555.48  —1555.55 
200 | —1576.36 | —1573.67 —1573.71 | —1575.52  —1575.61 | —1574.12  —1574.17 | —1572.54 —1572.61 
202 | —1592.19 | —1590.04 —1590.06 | —1591.58  —1591.65 | —1590.48  —1590.53 | —1589.35 | —1589.40 
204 | —1607.51 | —1606.03 —1606.05 | —1607.15  —1607.20 | —1606.50 —1606.53 | —1605.87 | —1605.91 
206 | —1622.32 | —1621.54 —1621.55 | —1622.11  —1622.13 | —1622.04  —1622.06 | —1622.00  —1622.02 
208 | —1636.43 | —1636.29 —1636.29 | —1636.00 | —1636.01 | —1636.80 —1636.81 | —1637.39  —1637.40 
210 | —1645.55 | —1644.26 —1644.23 | —1643.54 —1643.55 | —1644.79  —1644.77 | —1644.19  —1644.20 
212 | —1654.52 | —1652.13 —1652.04 | —1650.98 | —1650.99 | —1652.65  —1652.59 | —1651.00  —1651.01 
214 | —1663.29 | —1659.88  —1659.74 | —1658.31  —1658.30 | —1660.35  —1660.26 | —1657.82  —1657.83 


TABLE IV. Binding energy Ep (MeV), quadruple deforma- 
tions (8, Bn, Bp), matter and charge radii r and ren (fm) 
of ??Rn calculated by PKO1 (upper panel) and DD-ME2 
(lower panel) with respect to max, in which Km,,,, = 3, 
Ap = 0,2,4,6,8, and the initial deformation Bo = 0.2. Notice 
that EZP = —1697.796 MeV [115] and &**? = 0.1269 [117]. 


Mmax Ep T Teh B Bn Bp 
13/2 | —1696.27 5.753 5.659 | 0.1241 0.1275 0.1189 
15/2 | —1696.89 5.755 5.661 | 0.1297 0.1330 0.1246 
17/2 | —1697.20 5.757 5.663 | 0.1334 0.1369 0.1280 
19/2 | —1697.28 5.758 5.663 | 0.1345 0.1381 0.1290 
13/2 | —1694.75 5.726 5.656 | 0.1237 0.1273 0.1182 
15/2 | —1695.17 5.728 5.658 | 0.1271 0.1305 0.1217 
17/2 | —1695.34 5.729 5.659 | 0.1287 0.1322 0.1232 
19/2 | —1695.36 5.729 5.659 | 0.1288 0.1323 0.1233 


TABLE V. Deformations (8, Bn, Gp), binding energies Ep 
(MeV), and matter radii r (fm) and charge radii ren (fm) of 
Ne at local minima given by PKOi (i = 1,2,3) and DD- 
ME2. Notice that EZ? = —160.65 MeV [115]. 


Ep B Bn Bp r Tch 

140.74 | —0.798 —0.791 —0.805 | 3.054 3.168 

PKO1 153.56 | —0.121 —0.117 —0.124 | 2.745 2.873 
—159.43| 0.561 0.554 0.567] 2.820 2.943 

140.12 | —0.817 —0.810 —0.824 | 3.081 3.194 

PKO3 153.13 | —0.121 —0.118 —0.124 | 2.755 2.883 
—159.26| 0.580 0.573 0.586 |2.835 2.958 

PKO2 152.86 | —0.080 —0.077 —0.084 | 2.726 2.855 
—155.56| 0.520 0.514 0.527 2.791 2.915 
DD-ME2 151.53 | —0.180 —0.176 —0.184| 2.764 2.895 
—156.26| 0.541 0.533 0.548} 2.813 2.940 


method [104] than the BCS method [118]. The local 
minima extracted from Fig. 2 are summarized in Table 
V. Within the range 6 € (—1.0,1.0), there exist three 
local minima in the PKO1 and PKO3 results: a stable 
and strongly deformed prolate, a soft oblate, and a high- 
lying and largely deformed oblate. For the formal two 
prolate and oblate shapes, PKO2 and DD-ME?2 present 
rather similar deformations, while the largely deformed 
oblate is not supported. 


From Table V and Fig. 2, one can see that the selected 
effective Lagrangians give similar binding energies on the 
oblate minima [8 ~ (—0.1, —0.2)] of Ne, particularly 
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for the RHF ones PKOi. When approaching the prolate 


ss minima, the ground state of ?"Ne, the deviations between 
s» PKO2 and PKO1 (PKO3) become more and more 
seo notable. 
sı the 7-PV coupling, show much better agreement with the 


Eventually, PKO1 and PKOS3, which contain 


experimental data of the binding energy of ??Ne [115]. 


s3 In order to understand the deviations between PKO2 


and PKO1 (PKO3), we show the evolution of neutron 
(left panel) and proton (right panel) single particle 
energies with respect to the quadruple deformation 
8 € [0,0.6) in Fig. 3, where the notation v[m]z is 
introduced to denote the orbits. It can be seen that both 
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FIG. 3. (Color online) Evolution of neutron (left panel) and proton (right panel) single particle energies given by PKO2 (dotted 
lines) and PKO3 (solid lines) for ?"Ne with respect to the deformation £ € [0, 0.6), in which the filled and open circles denote 


the Fermi energies Er given by PKO3 and PKOJ2, respectively. 


neutron and proton single-particle energies show similar 
systematical evolutions with respect to 8. Following 
the shape evolution, both neutron and proton spherical 
shells N/Z = 8 are continuously reduced and eventually 
the deformed ones N/Z = 10 emerge at rather large 8 
values. In general, the orbits, which are branched from 
the same degenerated spherical orbits nlj (j > 3/2), will 
deviate from each other with respect to the deformation 
B. As seen from the PKO2 results (dotted lines), the 
orbits 2[1/2]+ and 1[1/2]— with the smallest m values, 
respectively branched from the spherical 1d5/5 and 1p3/2 
states, tend to be more and more bound, whereas the 
ones 1[5/2]+ and 1[3/2]— with the largest m values 
become less and less bound, along the shape evolution 
from the spherical to prolate. This in fact reflects the 
typical deformation effects. 


However, comparing to the results given by PKO3 
which contains the z-PV coupling [78], the situation 
can be notably different. As shown in Fig. 3, PKO3 
(solid lines) and PKO2 (dotted lines) present different 
systematics on the shape evolution of the valence orbits 
branched from the degenerated spherical one 1d5;5. In 
order to understand the effects of the 7-PV coupling, Fig. 
4 shows the splittings between valence neutron orbits 
2[1/2]-- and 1[3/2]+, namely AE = FEAj2)4. — E»n 2) 
(MeV), with respect to the quadruple deformation £, 
and the inset shows the effective neutron pairing gap A; 
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FIG. 4. (Color online) Splittings (MeV) between neutron 


valence orbits 2[1/2]4- and 1[3/2]+ given by PKO3 (solid line) 
and PKO2 (dotted line), namely AE = Fy/3/24 — Eoi1/2)+, 
with respect to the deformation 6 € [0,0.6) for 20Ne, where 
the grid pattern denotes the contributions from the «-PV 
coupling in PKO3. The inset shows the neutron pairing gaps 
An (MeV) as functions of 8 given by PKO3 and PKO2. 


(MeV) as functions of 8. It can be seen that following the «» PKO2 (dotted line), and the enhancements from PKO2 
shape evolution from the spherical to the prolate, PKO3 «0 to PKO3 are almost fully due to the contributions from 
ei (solid line) presents notably larger splittings AF than 6.2 the 7-PV coupling (in grid pattern). 
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FIG. 5. (Color online) Proportions (in percentage) of the 


main components in expanding the neutron 2[1/2]4- orbit for 
20Ne with respect to the deformation 8 € [0, 0.6], namely the 
1ds;5, 281/2, and 1da3;5 components. The results are given by 
PKO3 (solid lines) and PKO2 (dotted lines). 


Qualitatively, one can understand such enhancement 
from the nature of the tensor force carried by the 
7-PV coupling in PKO3 [78, 90], which leads to 
repulsive/attractive couplings between the « components 
with the same/opposite signs [119]. As a supplement, 
we show in Fig. 5 the proportions (in percentage) of the 
main expansion components [see Eq. (29)] as functions of 
the deformation 8 for the neutron orbit 2[1/2]+, namely 
the & components 1d5j2 (& = —3), ldajy9 (k = 2), 
and 25,/5 (k = —1). For both PKO2 (dotted lines) 
and PKOS3 (solid lines), the negative & component 1d5/2 
dominates the expansion of the Dirac spinor (more than 
60%). In fact, for another two neutron orbits 1[3/2]+ and 
1[5/2]--, which are not shown for simplicity, the negative 
k component 1d5;5 plays a compressively dominant role 
(more than 9096). Similar situation can be also found 
in expanding the proton orbits. If looking at the core 
configurations of both neutron and proton (N,Z = 8), 
only one orbit 2[1/2]— can be dominated by the positive « 
components 1p1/2 (x = 1), and the other three (1[1/2]+, 
1[1/2]-, and 1[3/2]-) contain mainly the negative & 
components, which indicates the core exhibiting as the 
negative & components in total. Thus, the tensor force 
components carried by the a-PV coupling in PKO3 
presents additional repulsive interactions between the 
core and valance orbits 2[1/2]4-, 1[3/2]+, and 1[5/2]4-. 

However, for the orbit 2[1/2]+, the proportions of the 
positive & component 1d3/2 (x = 2) are notable, which 
are largely enhanced from PKO2 to PKO3 following 
the shape evolution. In contrast to the negative « 
component lds5/2, the tensor couplings between the 
positive & component 1d3/2 and the core, that behaves 
like negative x components, are attractive, which tends to 
make ??Ne more bound. Comparing to the orbits 1[3/2]+ 
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and 1[5/2]--, the repulsive tensor couplings between the 
orbit 2[1/2]+ and the core are then notably reduced 
by enhanced proportions of the positive & component 
1d3/2, from which one can understand well the enhanced 
splittings AE from PKO2 to PKO3 in Fig. 4, as well as 
more bound ground state given by PKO3 than PKO2 in 
Fig. 2. 

On the other hand, via the a-PV coupling, the 
core of ??Ne can be also polarized differently by the 
valence nucleons in the PKO2 and PKOS calculations, 
manifested as notably different compositions of the 
valence orbit 2[1/2]+ in Fig. 5. It is thus expected 
that the mean field described by PKO2 and PKO3 can 
be notably different. As the direct evidences, not only 
presenting different compositions of the valence orbit 
2|1/2]4-, PKO2 and PKO3 also give distinctively different 
systematics on the shape evolutions of the orbits 1[3/2]4- 
and 1[5/2]+, see Fig. 3. Coincidentally, one may 
also find some similarity on the shape evolutions of the 
orbits 1[3/2]+ and 1[5/2]4-, and those of the 2s,/2 and 
1d3/2 proportions for the orbit 2[1/2]--, see Figs. 3 and 
5. To some extent, such similarity might illustrate the 
consistent relation between the systematics of the orbits 
1[3/2]+ and 1[5/2]+ and the core polarization by the 
valence nucleons occupying the orbit 2[1/2]+. It shall be 
noticed that the shape evolution of the orbits 1[3/2]+ and 
1[5/2]+ cannot be attributed directly to the role of the 
tensor force, because both orbits are similarly dominated 
by the negative x component 1ds5/2. Moreover, the 
pairing correlations may not play a substantial role in 
interpreting the systematic difference between PKO2 and 
PKOS3, since both show similar pairing effects, as seen 
from the inset of Fig. 4. 


IV. SUMMARY 


In this paper, the axially deformed relativistic Hartree- 
Fock (RHF) model is established by utilizing the 
spherical Dirac Woods-Saxon (DWS) base that can 
probably describe the asymptotic behavior of wave 
functions. The formalism of the axially deformed RHF 
model based on spherical DWS base is presented in 
details, as well as the relevant space truncations. Taking 
the light ??Ne, midheavy ?9Fe, and heavy Pb isotopes 
as the candidates, the reliability of the expansions are 
illustrated, as well as the verification of the accuracy of 
the deformed RHF codes. 

The axialy deformed RHF model based on the 
spherical DWS base is applied preliminarily to study 
the structure properties of ??Ne, a typically deformed 
nucleus, with the RHF effective Lagrangian PKOz and 
RMF ones DD-ME2. It is found that the RHF 
Lagrangians PKO1 and PKOS3, which contain the m- 
pseudovector coupling, improve the description of the 
binding energy of ??Ne. Furthermore, the discussions 
on the systematics of the shape evolution of single- 
particle energies indicate that the 7-PV coupling, mainly 


7 the tensor force component, may present notable effects 
n2 in determining the shape evolution of nucleus, as well 


m3 as the polarization of the core of nucleus. 


Thus, as 


74 perspectives, it is valuable for investigating the effects of 
ns the Fock terms, especially the 7-PV and p-T couplings, 
1e in determining the structure properties of the deformed 
n; nuclei, which deserves further systematical investigations 
ns With the existing RHF Lagrangians PKOi and PKAI. 
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Appendix A: Energy functionals and self-energies in various coupling channels 
1. Compounded symbols 


In the energy functionals, it is more convenient to express the expansion of the coupling strengths as 


golpo) —vV2n V g% (r)Yi,o(9, 9); (A1) 
Xp 


where gg corresponds to go, gw, Jp, and fr. In practice, we performed the integrations with respect to the angle 
variables Q = (V, y) since the expansion of the wave functions Wyrm is carried on the spherical Dirac Woods-Saxon 
(DWS) base nm. Such integration contains the Harmonic functions given by the expansions of the propagators 
(30), the coupling strengths (A1), and the couplings between the spherical spinors, 


laaa 
Von | dQY,,, (Q)Y, u, (QY 0(9) 7g As OX Dias ie (A2) 


where u = pa = fy, and (Aa; ua), (^y, Hy) and A, denote the terms due to the couplings between the spinors, the 
expansions of the propagators and the coupling strengths, respectively. As an abbreviation, we introduced the symbol 
O to denote the above integration, 


Aye Los 6 caA5A0 A50 
d z(-1)" 7 AqAy Ap ROCA (A3) 


V2 


For the couplings between spherical spinors, we introduce the following symbols 2 (2) and Q (Q) as 


W ook & 

A > ^ Y—1~7LO A 

imd 72A (ga Un edet (A4) 
M Mi 

DAR Is e[- 5 M uk (A5) 
A ji tl,—4 mam A0 JM a JIM 

Cian =(—1) t3 V31 jalil SCA i0CXMo h lp ANUS ie (A6) 

J y 51 
LA. Mi 
Ol iiss EU B D (AT) 


In the symbols 2 and J, jj = ma — m, and ji = ma + mı, and for Q and Q, one can find u +o = m2 — m, and 
ii d- 0 — m» 4 m. With these symbols, the couplings of the spherical spinors can be simply expressed as 


—1)mı +3 = 
gige oa = gua auus Op 08 = ———— YONG Vogue, TAB) 


KIM] KAM K17mj;K2712 


Jim; J2m2 On Jimi j2m2 [Dr 
Xa 
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oi Where eg is the covariant spherical base vector with c = —1,0,+1. Notice that the quantity « denotes the combination 
o of the quantum numbers (7,1). 


983 2. Energy functionals and self-energies from the Hartree terms 


For o-S coupling and the time component of w-V one, the self-energies from the Hartree terms can be expressed as 


Ey (r j= ) a r) f rar 2 Ora R$ a (rr) ) exert r^p" (r^), (A9a) 
XM 


AL 


xo tr HEN nf fdr Ds Onin FE 9 (rr DDAN gt r')p, 2 (r^). (A9b) 


MAG 


Thus, the relevant energy functionals read as 


2 
EP 24 = rdr X` p (r)E3 (r), po r?dr » P$ COPIA (r). (A10a) 
Xa 


; : . A, . A, 
oa For the time component of p-V coupling, such expressions can be obtained by replacing gà» (g^) with gp? (gp”) and 
985 pS (pp) with 2 (pz), where py = prs — eB For the Coulomb field (A-V coupling), these expressions can be 
oe deduced by setting Àj = à, = 0 and replacing the nucleon density p» as the proton one pp. For the Hartree terms of 
o; the space components of the vector couplings, as well as the 7-PV couplings, the contributions are derived as zero. 


Since the coupling strengths go, gw, and gp are density-dependent, the variations of the energy functionals (A10) 
may lead to the additional rearrangement terms as 


Ap , 
D,A A0 go A 2 A0 Ap AL 
ED (p =I > V) ide = ey] fr dr R$ y (nr) XO e». les? (A11) 
Ay Apra Po ALAL 
yo TO x Ab X 
D;,A 0 w A 0 i 
DPM (n) = 4.2 5 D> S, F z ay]. f rar! $n) DO, [e an]. (A12) 
— Ay ApÀa Po MAL 
988 3. Energy functionals and self-energies from the Fock terms 


For the Fock terms, the expressions are much more complicated than the Hartree terms. Notice that in deriving 
both the Hartree and Fock terms, we have set the angular momentum projection m to be positive for convenience. 
Thus, we need to revise the expansion of the spinor Yyrm. In principle, one can set the time conjugation partner as 
Vys—m = Prams P, being the time reversal operator, and the expansion of Yyr-m reads as, 


Purr -— qve ree (A13) 


o» Where, by definition, ,;4, and Yyk-m share the radial components G;, and Fik. For the Hartree terms, it does not 
oo bring additional complexity. While for the Fock terms, one has to carefully treat the couplings between the orbits m 
oot and —m/ (m, m’ > 0). It can be found that the partner (—m, —m/) present identical contributions to the one (m, m"), 
»» and the partners (m, =m’) and (—m,m’) provide identical contributions. 


'To express the contributions of the Fock terms in compact form, we introduce the symbol D for o-S coupling and 
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the time component of the vector ones as 


ad pi, A m 1 F A A, P — 
prow? s aH aH ^un ^ad aH Ayh 
KiKam;k. k5m/ ED [um d IER AU AL st: a Ki d uu AU A 3 , (A14) 
adh 
A ABA Ayit— 1 px M -X 
prey Enos aH ab Ay Q^ 4 Grab al^ ^y 
kikom;k kam — | kimrim Z Kom am OAD » On, ^b 2 Ki Lat V komikan ASA KR A! sb (A15) 
EU 
ABA Ay; —+ 1 M A ALB 
prey ZS aH ab Ayb 4 Grab aH Ayi 
kikom;k kgm — 2 [P kimiri m Arenlor AA B ^b cm —k'm' DEO oF At ar (A16) 
dad 
AApAD Ayia 1 | ogan Aun Ayh Aaf 2^ Ay 
RiKam;k. kgm! T9 2 Kim; nimi 2 Km; PS R 555 Ab + 2 Kim; "T Kym’; mU RA Ap]. (A17) 
rar, c 


Notice that the terms contains Z symbols correspond to the contributions from the partners (m, m^) and (—m, =m’), 
and those containing the symbols 2 for the partners (m,—m’) and (—m, m^). Thus, the nonlocal self-energies of the 
o-S coupling can be expressed as 


: 1 At arp, Ayitt 
K1,K2;0 __ ++ / A P [ml o pp? 
b pula =} on X Ore! X Riiki atm (r, r ) X Ig” (1) go (r ) ) Ryd, T,T Poo !:k1Ko m? (A18a) 
nm! KA KR, prt 
. 1 X AAÀpALAyit— 
K1,K2;:0 __ F— / A p [nl c / Dp^p Ms 
Noam = 22 X ) rr ; Rit rh, n'm’ (r, T ) X Ig” (r)go (r ) ; Ryd, TT ee nina (A18b) 
nm! KA KR, Apr A 
f 1 rv a rpd, Ay; — + 
K1,K2;0 __ —+ / A p Ím! c / PAprAys 
Ae am Ics on X Orr! ` Ri kin (r, r ) X go" (r)ge (r ) X RB T.T IP v tuti (A18c) 
mm! KA KR, ApAp A 
: 1 X a rpX, Ayi —— 
K1,K2;0 __ == 1 A Dp (Vl c / PprOy? 
A pas =+ on E Ore! ; Rairi ami (r,r ) X Ig” (r)go (r ) X Ryd, TT LC icm (A18d) 
mm! KR, Aprp A 


where the sum over v’ has been absorbed into the nonlocal densities R, and the factor 6,, indicates that the orbits 
mm and m'm’ should be both neutron or proton ones. In terms of the nonlocal self-energies, the energy functional of 
the o-S coupling reads as, 


K1,K2;7 K1,K2;0 
drd Yarm bg ans Gira Ai 
rar 2,5 5 ( Giri Fini), XELLE KX kiko X ’ ( 9) 
K i K2 Gyrm Firm r,r! 12 J p 


where v? (€ [0,2]) is the occupation number of the orbit (vnm). For the rearrangement terms, the contribution to 
the self-energy by can be similarly expressed as 


K1,KR2;0 K1,K2;0 
x5 Ad dpi (Gi FE mou PE M Gina A20 
T 2 po Lust Fini), C) ne Qe E ? ( ) 
K1K2 G, nm, Aa Fyrm,Xa r,r’ tha} p 


where the terms P and Q read as 


K1,K2;0 1 + 1 age” (r) ot c DS NF 
Pom mn 5 Orr! 5 RE ki nm (r,r P» UM (r i FS AGT LAE rm (A21a) 
nm! KA RS »^5 y 
K1,K2;0 1 H / g>” (r) uF pe Ayit— 
Erm Ad => on 2. Orr! 5 Rd id nm! (r, T p» Op" (r) me n2, RX, Ay (r, T ND Mum 'ik1Kom? (A21b) 
KRG prp 
K1,K2;0 1 1 gi" (r) (r) AApàp Ayit 
Ce cm. da —= on 2. Orr! 5 RR am! d T p» 8 Op" (r) me n2, RX, Ay (r, T ND R^ Km! "ium (A21c) 
+ KIKS p 1 
TE 1 o >P | T) X T ApBALAyi-— 
Qu A =+ on 5 Orr! 5 Rein, n'm’ (r, r’) ` ae, Jo" 02, Ri, Ay (r, r ND ki kam ae (A21d) 


mm! Kh Kh, Apà5 9p," (r) 
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ə3 For the time components of the vector couplings, similar expressions can be obtained by replacing the expansion 
s4 terms of the propagator RS a, and the coupling strengths ga" in Eqs. (A18), (A21) with the corresponding ones, 
os and additionally the plus signs in the Ya, Xr, Pg and Qr terms should be changed as the minus one. It should be 
os also noticed that for the isovector p-V coupling, the factor ô+ shall be replaced by (2 — 6,,/). For the Fock terms 
ə7 of the Coulomb interaction, there is no arrangement term and for the non-local self-energies [see Eqs. (A18)] only 
98 Àp = À, = 0 terms remain. 

For the space component of the vector couplings, we take the w-V coupling as examples, and the others can be 

similarly deduced. To write the expressions in a compact form, we introduce the symbols B as 


Apr, Ayit 1 [Xa Auc ^ae rN fio 
Bis simis iom = 5 yee ee e T x, + OEY TRO ME SE Mp (A22) 
aro | 
BrpAp Ayit— 1 [ nano No año ANI Bo Ayü | 
ee —9 5 OO au En Ab + QE m! En NE A VOX Ap | ’ (A23) 
Aare | 
QeApAy-t 1 » [oan Agno MH o Aaüo Aao ] 
Po seh cae nikai ~9 E n me a Pa B X X, F QO, /^m!,—&1m Qu RAP. SEU A, | ’ (A24) 
AaMj 
G^pApiAyi—— . 1 [| nano No WH O año Ni po Ay | 
P estet sm ise T9 > Oi mirim nena Ags Ha X, nd oe m/ "xu. cui AU ALI (A25) 
AaAMgm 
Thus, the nonlocal self-energies can be expressed as 
, 
: 1 x ANN Ayi t+ 
Ki1,K2;0 —— X y Eon d » A P (pl X ) Ww 1 p yi 
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; i X APA Ayit— 
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nm! KL RS Apr, À 
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Ki,K2;0 __ E RIS 1 A P (al wW 1 PAprAYys 
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Here, we use the bold type w to represent the space component. In terms of the above nonlocal self-energies, the 
energy functional from the space component of the w-V coupling reads as 


K1,K2;02 K1,K2;02 
Yanm Yanm Gira A 
drdr' Lm 2.4 (Gir Fins), KAR yi Raw F i ( 30) 
"m Gyrm Fyrm rr! iR2 y 
For the rearrangement terms, the contribution to the self-energy can be similarly expressed as 
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where the P and Q terms read as 
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o» For the other vector couplings, the expressions can be obtained by replacing the expansion terms of the propagator 
1000 and coupling strengths with the relevant ones. In addition, for the isovector p-V coupling, one needs to replace the 
11 isospin factor 0. by (2 —6,,/), and for the Coulomb interaction only the A, = n — 0 terms remain. 
For the 7-PV coupling, the situation becomes even more complicated. First there exist gradient operations over 
the propagator, which can be expressed as 
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where the radial terms read as 
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To have compact expressions, we introduce the symbols P and P as 
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where pto=p'+o' 2 m — m and i-o — ji +o’ — m 4 m/. With the defined symbols, the self-energies from the 
m-PV coupling can be expressed as 
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In terms of the non-local self-energies, the energy functional from the 7-PV coupling reads as 
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Similarly the rearrangement term X5 ^^ can be expressed as 
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where the terms P and Q read as 
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Besides, the contact term is introduced to compensate the zero-range term in Vj" " [see Eq. (A37)]. The Hartree 


contributions from the contact term are derived as zero, and the Fock contributions read as 
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To express the contact term compactly, we introduce the symbols %2 as 
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Thus, the energy functional of the contact term can be expressed as 
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where the self-energies read as 
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The rearrangement term in the self-energy can be derived as 
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with the P and Q terms as 
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